High-pressure phase transitions in BiFeOs: hydrostatic vs. non- hydrostatic conditions 
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We report high-pressure x-ray diffraction experiments on BiFeOa (BFO) single crystals in 
diamond-anvil cells up to 14 GPa. Two data sets are compared, one in hydrostatic conditions, 
with helium used as pressure-transmitting medium, and the other in non-hydrostatic conditions, 
with silicon oil as pressure-transmitting medium. It is shown that the crystal undergoes different 
phase transitions in the two cases, highlighting the high sensitivity of BFO to non-hydrostatic stress. 
Consequences for the interpretation of high-pressure structural studies are discussed. 
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I. INTRODUCTION 

Bismuth ferrite (BiFeOs - BFO) exhibits at ambient conditions both ferroelectric and magnetic order and has 
become the archetype of an intrinsic multiferroic, studied for both technology-driven and fondamental research 
From a structural point of view, BFO crystallizes at ambient conditions in the perovskite structure with space group 
i?3c. This structure can be derived from the cubic parent perovskite by rotations (tilts) of the oxygen octaedra around 
the [lll]c direction relative to the parent cubic cell {a~a~a~ in Glazer's notations) and displacements of the Bi'^+ 
and Fe^+ cations along the same [lll]c direction. 

Here we are interested in the effect of high-pressure, a parameter which has been much less explored in the past, 
when compared to the wealth of temperature-dependent measurements, and this although Samara's pioneer work 
0-111 amplyillustrated the usefulness of this parameter. More recently, the parameter high-pressure has attracted 
more interest @, as exemplified by studies on classic ferroelectrics ferroelastics [lol - [T3| . piezoelectrics 

[iTj etc. Although these recent studies have led to a better understanding of the effect of high-pressure on cation 
displacements in ferroelectrics or octahedra tilts in ferroelastics, there is only little work on materials that present 
simultaneously those two competing structural instabilities. BFO is one of the very few perovskites that exhibits both 
cation displacements and octaedra tilts at ambient conditions, making it an ideal candidate for this field of study. 
One of the consequences of this complex competition in BFO is a rich pressure-temperature phase diagram which 
is still not completely understood ; not only the high-temperature regime but also the effects of a high hydrostatic 
pressure remain controversial [l[ . In this study, we will focus on the effect of low and medium pressures (up to 14 
GPa) at room temperature, aiming at a better understanding of the literature discrepancies, which are summarized 
in the following. 

Powder x-ray diffraction under hi gh-p ressure has been performed over the past few years by several groups with 
various outcomes. Gavriliuk et al. |l8l . [l9j and Zhu et al. [l^ have reported no phase transition in this pressure 
range. Haumont et al. |2l| have identified two phase transitions at 3.5 and 10 GPa with the phase sequence i?3c — > 
C2/m — > Pnma. Belik et al. [l^ have later confirmed the phase transition at 4 GPa and identified an additional 
transition at 7 GPa. They have identified the two intermediate phases as orthorhombic instead of the monoclinic 
C2/m, with a transition sequence i?3c — > Ortho I — > Ortho II — > Pnma, and also pointed out reversibility issues. 
In parallel, Raman spectroscopy has also been used to investigate phase transitions [23l. [2^. Both studies revealed 
two phase transitions in the low-pressure range, but at somewhat different pressures, about 3 and 9 GPa. 

This variety of results reported in the literature suggests that the structure of BFO with its competing instabilities 
is particularly sensitive to the quality of the sample and experimental conditions (stoichiomctry, homogeneity and 
isotropy of the stress field etc.). In order to evaluate the effect of non-hydrostatic pressure conditions as one source 
of the controversy, we compare in this work two high-pressure single-crystal x-ray diffraction experiments carried out 
in diamond anvil cells with different pressure-transmitting media. This comparison suggests that the quality of the 
hydrostatic conditions is one source explaining the variety of results in the literature on BFO. 
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II. EXPERIMENTAL ASPECTS 

BFO single crystals were grown by the flux method as described elsewhere [2^ . The samples used for the experiments 
were polished to a thickness of about 10 |J.m with a lateral extension of 10 to 30 [a.m. No domain structure was visible 
before the experiments, neither under polarized light nor in the initial diffraction patterns. All experiments were 
performed in diamond-anvil cells (DAC). The diamonds have the Boehlcr-Almax design with a cullet of 600 |xm. The 
pressure chamber was sealed by a stainless steel gasket pre-indentcd to a thickness of about 60 |xm. 

X-ray diffraction experiments were performed on the ID27 and ID09A beamlines at the ESRF. At ID27, the beam 
was monochromatic with a wavelength of 0.3738 A selected by an iodine K-edge filter and focused to a beam size 
of about 3 |a.m. At ID09A, the beam size is about 20 [a.m and the wavelength (0.4144 A) was determined from the 
calibration using a standard silicon powder. The signal was collected in the rotating crystal geometry on a CCD 
detector with —30° < w < 30° or —24° < w < 24°, depending on the opening of the cell used. A precise calibration 
of the detector parameters was performed with a reference silicon powder. The diffraction patterns were indexed 
with a home-made program based on the Fit2D software [26| . The refinement of the lattice constants from the peak 
positions was performed with the program UnitCell (27l |. 



Loading axis 
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FIG. 1: Sketch of the diamond-anvil cell (DAC) with the orientation of the sample relative to the loading axis of the cell. The 
sample is immersed into a pressure-transmitting medium (PTM): helium and silicon oil for hydrostatic and non-hydrostatic 
conditions, respectively. 



For a measurement in best hydrostatic conditions, helium was used as a pressure transmitting medium (PTM). 
Helium remains liquid up to 12.7 GPa. As such, it provides perfect hydrostatic conditions below this limit and remains 
very good up to very high pressures [28l [29| . For the measurement in non- hydrostatic conditions, we used silicon oil 
as a PTM, which is a very viscous fluid that freezes rapidly. Its hydrostatic limit has been reported to lie around 3 
GPa, as determined from the analysis of the ruby luminescence lines [2^, or even 0.9 GPa from the measurement of 
diffraction line broadening induced in a quartz reference crystal [2^. Above this limit, the stress can no longer be 
regarded as hydrostatic and a so-called deviatoric component adds to the hydrostatic stress field. With the z-axis 
along the loading axis of the cell, the total stress field a can be written [s^l : 




where P and ud measure the magnitude of the hydrostatic pressure and deviatoric stress respectively. The deviatoric 
stress has a rotational symmetry around the loading axis of the cell, with a stronger compressive stress along the 
axis. In an experiment with silicon oil as PTM, Zhao et al. (s^l have estimated this stress to be of 1 GPa at 10 GPa 
approximately. Although our setup is not strictly identical, this value can be considered as a reasonable estimate for 
the maximum deviatoric stress in our experiment. 

For a non-hydrostatic stress field, the effect of high pressure is expected to depend on the orientation of the crystal 
with respect to the loading axis of the DAC. In that case, the loading axis is directed perpendicular to a natural face 
of the as-grown crystal, close to [101]c. The setup with the orientation of the crystal is sketched in figure [TJ 

In all the experiments, the fluorescence of ruby was used as a pressure gauge [31|. Note that as the non-hydrostatic 
stress sets on, the pressure measurement from the ruby line becomes slightly inaccurate, but it remains within negligible 
uncertainties [30| and does not alter the conclusions presented in this study. 
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III. RESULTS 
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FIG. 2: Diffraction patterns observed for the different phases. For clarity, all patterns are taken from the same experiment, 
namely with non- hydrostatic conditions. Patterns (a), (b) and (d) are qualitatively identical when the conditions are hydrostatic. 
Pattern (c) was only observed in non-hydrostatic conditions. The indices given in pattern (a) are relative to the parent cubic 
cell. 



We start with the description of the phase transitions observed for hydrostatic pressure conditions. At ambient 
pressure, the crystal is in a single domain state in the familiar rhomboedral structure with space group i?3c. The 
diffraction pattern integrated over the full angular range is identical to the pattern shown in figure [2] (a). Between 
room pressure and 4 GPa, the rhomboedral distortion, as measured by the rhomboedral angle a, descreases gradually. 
At 4 GPa, the crystal undergoes a first transition to a new phase which we label phase 11. In this phase, the crystal 
exhibits a more complex diffraction pattern, with the emergence of new superstructure reflections (figure [5] (b)). In 
addition, the 3 to 4-folding of many reflections reveals the presence of a complex domain structure. The pattern can 
be indexed within the monoclinic C2/m cell proposed by Haumont et a/.[2ll|. This cell has a comparatively large 
volume {Z = 12), with lattice constants = 17.602(1) A, 5^ = 7.762(1) A, = 5.510(1) A and /3 = 108.195(1)° at 
4.3 GPa. The unit cell vectors are sketched in figure [3] (a) . For this cell, the pseudo-cubic lattice constants plotted 
in figured] were calculated as Opc = am/2-\/5, fopc = &m/2 and Cp, = Cm/v^- At 7 GPa, a phase transition to a new 
high-pressure phase (phase III) is observed. The transition is evidenced by the evolution of the peak positions but is 
not associated with the emergence or disappearance of peaks in the diffraction pattern. The pattern could therefore 
be indexed within the same unit cell, but with sudden changes of the Opc and the Cpc lattice constants respectively 
(figure E]), while the volume remains almost constant. Finally, at 11 GPa, the crystal undergoes a transition to the 
Pnma orthorhombic structure, with a diffraction pattern shown in figure [2] (d). The lattice vectors are recalled in 
figure El (c) . 

The phase sequence then reads R3c — > Phase II — > Phase III — > Pnma. The transition pressures are consistent 
with the results by Belik et al. [l^l and Haumont et al. [2l[, except that the more subtle transition Phase II — > 
Phase III, only associated with the jumps of the Op, and Cp, lattice parameters, was overlooked in (2ll |. In addition, 
the lattice constants in the Pnma phase match very well the values reported in ref. (2l| from powder x-ray diffraction 
carried out with hydrogen as PTM, which confirms the hydrostaticity of the stress field. The precise identification of 
space groups for phases II and III is beyond the scope of this paper, but we recall that Belik et al. have proposed two 
orthorhombic phases for phases II and III, also sketched in figure [3] (b). This proposition is similar to the monoclinic 
C2/m adopted here for the calculation of lattice constants in the sense that both cells have the same volume, but 
differ by the presence of the monoclinic distortion. 
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We now turn to the observations made in non-hydrostatic pressure conditions. The evolution of the rhomboedral 
angle a in the two experiments at hydrostatic and non- hydrostatic conditions are compared in figure 2] (insert). The 
distortion decreases more slowly in the non-hydrostatic case, reflecting the onset of the deviatoric stress even at low 
pressures. The transition i?3c — s^Phase II is shifted to higher pressure by almost 1 GPa. Between 5 and 7 GPa, 
phase II remains stable and is identical to the phase observed under hydrostatic conditions. The transition to phase 
III occurs at a pressure that is very close to the 7 GPa observed in the hydrostatic case. However, at 8 GPa, the 
crystal exhibits a phase transition to a distinct phase, labelled Ill-b, that was not observed in hydrostatic conditions. 
The diffraction pattern is shown in figure[2](c). This pattern could be indexed in a monoclinic cell with Om = 9.560(2) 
A, 6m = 5.35(2) A, Cm 9.547(3) A and /? « 109.44(3)° at 10.2 GPa, with lattice vectors sketched in figure [3] (d). 
In this experiment, the short b axis lies almost in the direction of the loading axis of the cell and bears the larger 
compressive stress. The crystal exhibits again a domain structure that makes the determination of the space group 
difficult, but it is interesting to note that this cell is metrically very close to the unit cells observed in other bismuth 
based perovskitcs: BiCrOs and BiMnOs [i^,!!^. We also point out that this monoclinic cell has Om « Cm and is 
therefore very close to a face-centered orthorhombic cell. Finally, at 11 GPa, the crystal undergoes a transition to 
the orthorhombic Pnma structure. It is noticeable that the crystal is then almost in a single domain state, with its 
c-axis directed along the loading axis of the cell. The lattice constants in the Pnma phase are very different from the 
parameters obtained in hydrostatic conditions (figure 01), which is in agreement with the expected larger compressive 
stress along the loading axis. A comparison of the relative volume V/Vq under both conditions (figure U) shows that 
the volume measured under non-hydrostatic conditions is larger than in hydrostatic conditions. 



IV. DISCUSSION AND CONCLUDING REMARKS 



Because non-hydrostatic stresses are unavoidable at high enough pressures in diamond-anvil cells, the influence of 
deviatoric stresses on powder diffraction experiments in general has been given a lot of attention (see e.g. [33 - [36j 
and references therein). Non- hydrostatic conditions are known to cause broadening, shifting, and even splitting of 
the diffraction peaks, even in the absence of any phase transition. In addition, grain-to-grain stresses may have an 
effect on the structure itself, and lead to different experimental results. For that reasons, it is recommandcd to check 
powder data against single crystal data for the identification of phase transitions. 

In this paper, we have reported two different experiments on BFO single crystals, carried out in hydrostatic and 
non-hydrostatic conditions. We have shown that the non-hydrostatic stress causes a significant shift in the first 
transition pressure, and more importantly a modification of the transition sequence, with the emergence of a new 
monoclinic phase Ill-b. The emergence of a distinct phase is a remarkable effect and shows that a non-hydrostatic 
stress field is a particularly critical parameter for high-pressure experiments on BFO. It is expected to be a major 
concern when the PTM becomes solid. From a practical point of view, great care should be taken with the PTM and 
the experimental conditions in general. 
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FIG. 4: Evolution under pressure of the pseudo-cubic lattice constants (top), rhomboedral angle (inset) and relative volume 
(bottom). The full and open symbols correspond to hydrostatic and non- hydrostatic conditions respectively. The stability 
domains of the various phases represented as coloured stripes are given for hydrostatic conditions, while the phase sequence 
observed in non- hydrostatic conditions is indicated by dashed lines in the volume plot. 



With this in mind, we suggest that the different PTM used in literature may be the dominating cause of the different 
transition pressures reported from Raman spectroscopy in ref. [2^ (Argon as PTM, sohd above 1.3 GPa) and infrared 
spectroscopy in ref. [21| (solid Csl as PTM). We note t he g ood agreement between the results obtained in this work 
on single crystals and the results presented in ref. [2l|, |22| on powder samples, all experiments being performed in 
very good hydrostatic conditions. On the other hand, non-hydrostatic conditions cannot explain all the discrepancies 
in the literature. In ref. [l^, H^l, the low-pressure phase transitions have not been observed although the PTM used 
(methanol-ethanol mixture and helium respectively) should ensure very good hydrostatic conditions in this pressure 
range. Other sources of disagreement, such as sample stoichiometry, have to be considered. 
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